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O ■ We study the problem of folding of the regular triangular lattice in the presence of 

a quenched random bending rigidity ±K and a magnetic field h (conjugate to the local 
^ ! normal vectors to the triangles). The randomness in the bending energy can be understood 

■ as arising from a prior marking of the lattice with quenched creases on which folds are 
favored. We consider three types of quenched randomness: (1) a "physical" randomness 
where the creases arise from some prior random folding; (2) a Mattis-like randomness 
where creases are domain walls of some quenched spin system; (3) an Edwards- Anderson- 
like randomness where the bending energy is ±K at random independently on each bond. 
The corresponding (i^T, h) phase diagrams are determined in the hexagon approximation 
of the cluster variation method. Depending on the type of randomness, the system shows 
essentially different behaviors. 
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1. Introduction 

The statistical properties of polymerized membranes, or tethered surfaces have been 
widely discussed in the past few years 0-^. A polymerized membrane is the two- 
dimensional generalization of a linear polymer Its energy involves both an in-plane 
elastic (strain) contribution and an out-of-plane (bending) one. At low temperature, such 
a membrane with bending rigidity is asymptotically fiat and its radius of gyration Rc 
increases as the linear internal dimension L of the surface p|-|TO[| . As a function of temper- 
ature, the membrane without self-avoiding interaction (phantom membrane) undergoes a 
crumpling transition from the low temperature flat phase to a high temperature crumpled 



phase [Rq ~ vhiL) |T^-[T^]. The mechanism of the transition, in particular the stability of 



the flat phase, is rather subtle and relies on the coupling between in-plane and out-of-plane 
deformation modes [|5|. 

Very generally, 2-dimensional membranes can be discretized into triangulations, whose 
faces are endowed with natural Heisenberg spin variables, representing the direction of 
the local normal vector to the surface Polymerized membranes, which have a fixed 
connectivity, then translate into statistical spin systems on the regular triangular lattice. 
The corresponding spin system is however involved, because the resulting spins are not 
independent variables. The constraint of being normal vectors to a surface causes a long 
range interaction between the spins, which stabilizes the ordered flat phase, as opposed 
to the case of the usual 2-dimensional unconstrained Heisenberg model, always disordered 



The above mechanism clearly indicates the subtlety of the correspondence between 
geometrical objects and spin systems, especially in two dimensions. In order to understand 
this, several simple models have been proposed. The simplest one is probably a square 
lattice model introduced by David and Guitter 0. The model is a discrete rigid bond 
square lattice, which is allowed to fold onto itself along its bonds in a two-dimensional 
embedding space. The constraint there makes the folds propagate along straight lines. In 
the presence of bending rigidity, the system is always flat. 

A triangular lattice version of this 2D folding model was then introduced by Kantor 



and Jaric [IT^. Describing the (up or down) normal to each triangle by a spin variable 



a = ±1, the model Hamiltonian translates into that of the Ising model with however 
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constrains on the spin variables. This results in several new features, and a totally new 
phase diagram, as compared with the usual Ising model []T^|TB[. Finally, the folding of 
the triangular lattice embedded in a 3-dimensional discrete space has been formulated as 



a 96- vertex model [|T^,[T^]. 

Constraints do not appear only in the physical degrees of freedom, like local normal 
vectors. If the membrane has disorder, the disorder itself can also be constrained in some 
cases. For example, if one folds a piece of paper and makes a crease, this will generate a 
spontaneous curvature along the crease. If one now crumples the paper randomly by hand 
H,|T9| , the generated spontaneous curvature will be directly related to the configuration of 
the normal vectors of the resulting crumpled configuration. For the latter to be accepted 
as a physical configuration, the normal vectors should also obey some constraints. The 
induced disorder (in this case the induced random spontaneous curvature) should thus 
obey similar "physical" constraints. 

Here we study a simple model of folding with such a "physical" quenched randomness. 
As a model, we use the triangular lattice with quenched random bending rigidity and clarify 
the importance of the "physical" constraints on the disorder. 

The paper is organized as follows. In section 2, we first discuss the general folding 
problem of the triangular lattice and recall some known facts about it. We then describe 
the precise type of randomness which we consider on the lattice and which we write as a 
Mattis-like spin system with constraints. In section 3, we describe the cluster variation 



method that we shall apply to the study of the thermodynamics of the system [pT| - p3| . We 
first describe the procedure for a general disordered system. We then restrict ourselves to 
the hexagonal approximation in which the clusters are made of a maximum of six triangles. 
Next, we analyze the pure (without disorder) system again, as a particular limiting case 
with trivial disorder. The results for the fully disordered system are presented in section 
4. After giving a few results following from a reduced variable analysis, we present the 
complete phase diagram of the system. In section 5, we study other variants of the disorder. 
In particular, we discuss the specificity of the "physical" constraint on the random bending 
rigidity. We present for each type of disorder the corresponding phase diagram. Some 
concluding remarks are gathered in section 6. 



2. Folding with Quenched Random Bending Rigidity 



In this section, we first recall the rules of folding for the triangular lattice [0. We 
then introduce disorder in the problem in terms of a random bending rigidity. 
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2.1. Folding of the triangular lattice: pure case 

Let us consider a regular triangular lattice which can be folded onto itself along its 
bonds. We allow only for complete foldings which result in two-dimensional folded config- 
urations. Each bond thus serves as a hinge between its two neighboring triangles and is in 
either one of the two states: folded (with the two neighboring triangles face to face) or not 
(with the two neighboring triangles side by side). A folded state of the system is entirely 
determined by the list of its folded bonds. In this definition, the folding process may cause 
self-intersections and the model corresponds to a "phantom" membrane. Also this does 
not distinguish between the different ways of folding which result in the same folded state. 



(a) (b) (c) (d) (e) (f) 

(g) (h) (i) (j) (k) 

Fig. 1: The eleven local fold environments for a vertex. Folds are represented 
by thick lines. One of the two (opposite) possible spin configurations on the 
triangles is also indicated. 

One can easily see that, among the 2^ possible fold configurations for the six bonds 
surrounding a given vertex, only 11 states are allowed, corresponding to actual foldings 
of the surrounding hexagon [^. These configurations are displayed in fig 1. It can be 



checked that imposing everywhere one of these eleven local environments is sufficient to 
define the folding consistently throughout the lattice. The folding of the triangular lattice 
is thus simply expressed as an 11-vertex model on the lattice. Note that, even if the folding 
is defined locally, its nature is highly non-local. Since all vertices in Fig.l have an even 
number of elementary folds, these folds form folding lines without endpoints. Moreover 
all "folded" vertices (vertices (b)-(k) in Fig.l) have at least one fold on the left half of 
the hexagon and one on the right. Thus folds are forced to propagate through the entire 



lattice 14 



Folding can also be expressed with spin variables cr^ = ±1 living on the elementary 
triangles, which indicate whether the i-th triangle faces up or down in the folded state. 
The spin variable changes its sign between two neighboring triangles, if and only if their 
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common bond is folded, i.e. folding lines are domain walls of the spin system. One can 
think of the spin as the normal vector to the triangle. We depict the corresponding spin 
configurations on Fig.l. Note that there are two spin configurations for each folded state, 
due to the degeneracy under reversal of all spins. 

Clearly, the only allowed vertex environments are those with exactly 0, 3 or 6 sur- 
rounding up spins. In other words, for a spin configuration to correspond to a folded state, 
the six spins ai(i = 1, 2, ■ ■ ■ , 6) around any vertex v must satisfy the local constraint [P 



= ai = mod 3 , (2.1) 

i around v 

since E„ = 2 (number of up spins) — 6 is a multiple of 3 iff the number of up spins itself is 
a multiple of 3. Folding is thus expressed here as a constrained spin system. 

The statistical behavior of this system has been extensively studied, using a transfer 
matrix formalism [p!3|JT5|] , the correspondence with a solvable three-coloring model [0 



and a cluster variation method []T6[. Introducing a bending energy term —Jaiaj between 
nearest neighbors and a magnetic field term —Hai, the following model Hamiltonian was 
considered, 

Ti-ising = -j'^CTiCTj - H'^ai . (2.2) 

This is nothing but the Ising Hamiltonian, which is here however coupled to the local 
constraint ( |2.1| ). In the folding context, the magnetization simply measures the projected 
area of the lattice (i.e. the algebraic area of the domain enclosed by its boundary) and the 
magnetic field can thus be interpreted as a lateral tension term. For convenience we will 
use the reduced coupling and magnetic field 

K = J/ksT ; h = H/ksT . (2.3) 

The phase diagram of the system is shown in Fig. 2. In order to characterize each 
phase, two order parameters are introduced. One is the usual magnetization, 

M=^{{J2-^ + J2^^)), (2.4) 

* A V 

and the other is the staggered magnetization 

^«t = ^((E^^-E^O)' (2-5) 

* A V 
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Fig. 2: Phase diagram in the (K, h) plane. Three first order lines h = 
hc{K), —hc{K) {K < Kc) and h — {K > Kc) separate the three phases 
M = 0, ±1 and meet at the triple point {K^ 0). The dashed line represents 
the transition line between the disordered folded phase M = 0, Mst = and 
the compactly ordered folded phase M = 0, Mgt ^ 0. 

where the sum is performed separately on triangles pointing up and down, dividing the orig- 
inal triangular lattice into two inter-penetrating sublattices. Here Nt is the total number 
of triangles in the system. Three phases exist: a completely fiat phase (M = ±1, Mgt = 0), 
a disordered folded state (M = 0, Mgt = 0) and a compactly ordered folded state 



(M = 0, Mst 7^ 0) |]T5|,|T^. Note that the fiat phase has a maximal magnetization |M| = 1 
and is indeed frozen in the pure completely fiat state with all spins aligned. There is no 
fiat phase with intermediate (0 < |M| < 1) magnetization. Three first order transition 
lines h = hc{K), —hc{K) {K < Kc), and h = {K > Kc) separate the three phases 
M = 0, ±1, with a triple point at Kc ~ 0.1 (estimated from either the transfer matrix 
T5[1 or the cluster variation |jT^ approach). For negative K, the transition between the 



disordered folded phase and a compactly ordered (antiferromagnetic) folded phase with 
staggered order parameter Mgt 7^ is found to be continuous at h = []T^ . This transition 
is represented by the broken line of Fig. 2, which intersects the horizontal {h = 0) axis at 
K ~ —0.284. We see here that, when compared with the usual unconstrained Ising model, 
the phase diagram has been strongly modified. In particular, it is now asymmetric with 
respect to K with the usual continuous ferromagnetic transition replaced by an abrupt 
first order transition to a completely ordered phase. 
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2.2. Quenched random bending rigidity 

What kind of disorder should one introduce in the above model (Eq.( |2.2|) )? Since 
we are dealing with folding of a phantom object with a two-dimensional resulting folded 
state, we cannot distinguish between ±180 degree folds, thus we cannot introduce any 
spontaneous curvature term. Disorder will appear here in the form of a random bending 
rigidity Kij. The most general Hamiltonian with random nearest neighbor coupling Kij 
and a constant external field h is (we drop the ksT factor, thus using reduced coupling 
constants) 

'^random = ^ — -f^ijC^iCTj — ^ CTi ■ (2.6) 

There are several possible choices for Kij . One possibility is the Edwards- Anderson model 
for spin-glasses [^, where Kij = ±K independently on each bond. Such a model is 



interesting as it may have a spin glass phase. In such a phase, spins are randomly oriented 
(M = and Mgt = 0) , but the following order parameter q takes a nonzero value 



(2-7) 



Here the upper bar stands for the quenched average over the randomness. We shall discuss 
this type of disorder in section 5. 

As has been discussed in the introduction, we are interested in another type of "physi- 
cal" disorder where the random bending rigidity has been generated by a prior irreversible 
folding of the lattice [^|T^. We have in mind the picture of a crumpled piece of paper 
marked with irreversible creases. The effect of this irreversible crumpling is to impose the 
corresponding crumpled state as the new ground state of the system. No frustration will 
occur as long as one does not strain the paper. The model Hamiltonian should thus have 
that crumpled state as its ground state and this information should be included in Kij. 
One natural choice for Kij with a "random" ordered phase and no frustration is that of 



a Mattis Model |^^. In this model, the bending rigidities Kij are functions of a set of 
random face variables r^: 

K,j = Kr.Tj. (2.8) 

The Ti = ±1 are "frozen" according to a specified probability distribution Pt{ti,T2 ■ ■ ■ , tatJ 
reminiscent of the first irreversible crumpling process. The total Hamiltonian is then given 
by 

^Mattis = -K ^ TiTjCFiaj - h^ai . (2.9) 
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As is well-known, if there is no constraint on spin variables, the following gauge transfor- 
mation 

(j[ o,T, (2.10) 

makes the above model ( |2.9|) in zero external field /i = equivalent to the pure system 
||20[] . KX, K = 00, the state with Oi = Ti is recovered as the ground state. 

In our model, the spin variables are constrained by Eq. and the above gauge trans- 
formation cannot be performed since it does not preserve the folding constraint. Moreover, 
the ground state ai = Ti can be reached only if the r variables themselves obey the folding 
constraint: 

^ r, = mod 3 . (2.11) 

i around v 

This type of r-configuration is what we call a "physical" disorder. We shall restrict our- 
selves to this type of disorder in the following sections 3 and 4. We will return to other 
types of disorder (such as the Edwards- Anderson model) later in section 5. 

From the above discussion, we understand the physical origin of the local constraint 
on the r- variables. Then what probability distribution Priji, t~2 - ■• , t^J should we use for 
them? Since the r variables obey the same constraint as the pure a system described in 
the first part of this section, we can take advantage of the solution of this pure system 
by simply assuming that the r distribution is described by a particular (appropriately 
chosen) point in the disordered phase of the phase diagram of Fig. 2. A natural choice is 
the point at the origin of the phase diagram [K = h = 0) since it then does not involve any 
energy parameter for the disorder and treats as equiprobable all allowed configurations of 
"physical" disorder. In other words, we shall take for the distribution Pt(t"i, T2 ■ • ■ , tatJ 
the density of the pure constrained problem at K = h = 0. Due to the constraint, this 
density remains non-trivial. 



3. Cluster Variation Method for the Disordered System 

In this section we explain in detail the hexagon approximation of the cluster variation 
method (CVM) generalized to a random system [^T]-^. In subsection 3.1, we explain the 
method in the general case. In subsection 3.2, we apply the CVM to the pure system 
[Kij = K) as a limiting case with trivial disorder. As mentioned above, this is also 
instrumental to give an explicit form for pj- which is necessary to tackle the fully disordered 
case. We also discuss several symmetry breakings of the model. 
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3.1. The CVM and its hexagon approximation 

The CVM is a closed-form approximation based on the minimization of an approx- 
imated free energy density functional, which is obtained by a truncation of the cluster 
expansion of the full free energy density functional appearing in the exact variational 



formulation of the problem [p5l ,|2 



Consider our spin system with A^^ sites and Hamiltonian (|2.6|)[p3| . The configura- 
tion of the random bond couplings Kij is specified by a probability distribution P{{Kij}). 
In our case of "physical" disorder where Kij = KriTj, we shall have: 

Pi{K,j = Kt.Tj}) = prin, ■ ■ ■ , TivJ. (3.1) 

The discussion below is however more general. 

In terms of a density matrix p(ai, a2, • ■ ■ , crTvJI-ftTi-,}) for each configuration {Kij}, 
we define the variational free energy associated with the Hamiltonian Tirandomlc, {^ij}) 
(from now on, we use the notation a = {cxi}) as 

Ti{K,j}) = [y]p(a|{K,,})[Kandom(a,{K.,}) + lnp(a|{K,,})]l . (3.2) 

L ' J min 

a 

The subscript min means that the above expression must be taken at its minimum with 
respect to p{a\{Kij}). This is the well-known variational principle. The minimization is 
performed at fixed {Kij} with the normalization constraint: 

J]p(cT|{i^,,}) = l. (3.3) 

The quenched free energy JF is then given by 

^=5^ P{{K,,})J^{{K,,}), (3.4) 

where the sum extends over all possible realizations of the disorder. Upon introducing the 
generalized density matrix 

p(a, {K,,}) = P{{K,,})p{a\{K,,}), (3.5) 

one can easily show that 

y. p(^,{^^.})[^random(a,{K,,})+lnp(a,{K,,})]l + ^Dis (3.6) 
L ' J min 



where the minimization is now on a density p(o", {Kij}) for both a and {Kij} with the 
constraint 



(3.7) 



The quantity Sr,is is a constant term depending only on the probabhty distribution for the 
disorder and reads 



(3i 



In such a scheme it can be shown that the quenched average of the expectation value 
< A(cr, {Kij}) > of an operator A{a, {Kij}) is given by 



< A(a, {K,,}) >= Yl < ^(^' {K^j}) > 

= Yl Pmin{(jAKrj})A{a,{K,j}). 



(3.9) 



where Pmin(o"5 {^ij}) is the density at the minimum of the quenched free energy functional 

m- 

The CVM is obtained by taking the thermodynamic limit Nt ^ oo and truncating the 
cumulant expansion for the entropy S = —J2a {Kij} /^(^' i^ij}) ^^Pi'^^ {^ij}) appearing 
in ( |3.6| ) to a set of "maximal preserved clusters" r^,z = l,2,---,r (and all their translated 
images). The variational principle will then be applied to the reduced density matrix 
Pri{<^^ {Kij}) associated with the maximal preserved clusters F^, i.e. the minimization will 
be performed on this reduced set of densities. 




Fig. 3: Labeling of the spins on an elementary hexagon. Each site i (= 
1, ■ ■ ■ , 6) supports a spin variable cr^ and a disorder variable r^. The reduction 
process from pe to p2 and to Pia{b) is also indicated. 
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In the hexagon approximation for the triangular lattice, the largest clusters appearing 
in the expansion are hexagons. Hereafter we restrict our presentation to the case of the 
Mattis like coupling Kij = KriTj. For the other cases, the generalization is straightforward. 
We introduce the reduced density matrix for a hexagon. 

PeicTi, era, 0-4, a5, CT6, ^i-, ^s, T4, rs, re). (3.10) 

The spins in the argument of pe follow each other counterclockwise in the hexagon, and 
the first one is on the A sublattice (see Fig. 3) i.e. is pointing up. This reduced density 
matrix represents the probability for one hexagon to have fixed values of a and r. It is 
normalized according to 

E/^6(M,{t}) = p.,6({t}), (3.11) 

{^} 

where is the 6-point probability for the disorder variable on a hexagon as obtained from 
the corresponding partial trace of Pt-(ti, • • • , tatJ in the thermodynamic limit. We assume 
that this 6-point reduced density is the same for each hexagon, i.e. that the distribution 
of disorder is translationaly invariant. 

We also introduce the site and pair density matrices Pia(b)(ci, ^1), P2(<7i, a"2, ti, T2), 
which are defined as symmetrized partial traces of p% by 

P2(cri,Cr2,Ti,T2) = - E [P6(<7l,Cr2,Cr3,Cr4,Cr5,Cr6,n,T-2,T3,T4,T5,T6) 

+ P6((7'3, (y-l, (y\, 0-4, 0-5, (76, T3, T2, Ti, T^, Tr,, Te) 
+ P6(cr3, Cr4, CTl, Cr2, C^S, C^e, T3, T4, Ti, T2, T5, Te) 
+ P6((''3,(''4,(''55(''2,Cri,0'6,T3,T4,T5,T2,Ti,T6) 

(3.12) 

+ P6(o'3, o'4, o's, o'e, o"!, (72, T3, T4, r5, re, Ti, T2) 
+ Pe(cri, cr3, cr4, cr5, CTe, cr2, T-3, r4, r5, re, r2)], 
Pi^(cri,ri) = P2((^l,cr2,Ti,r2), 

PlB(cr2,T2) = E P2((7'l,Cr2,Ti,r2). 

0"l,Tl 

Here we have introduced two site density matrices, p\A and pis, corresponding to the two 
inter-penetrating sublattices in which the triangular lattice can be divided (see Fig. 3). 
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After the appropriate truncation of the cumulant expansion of S at the level of hexago- 
nal clusters, we get the approximate CVM quenched free energy per hexagon as a functional 
of peiWi}, {rj) only (by imphcit use of Eq.( |3.12| ))[p5|, 



f{p6{{(Ti},{Ti})) = -3i^Tro-,^(rir2cricr2P2(<7z,a2,ri,r2)) 
- /iTr^(cTipiA(o-i, Ti) - hTr^{a2PiB{cr2, r2) 

+ Tro-,^(p6lnp6) - 3Tr(^,r(p2 lnp2) + Tt^^ripiA^T^ Pia) + Tro-,^(piB Inpis) (3.13) 



+ Tr, 



A,({rJ)(TY,P6({^^}, {r^}) - PrA{r^})) 



to be minimized with respect to P6({o"i}7 Here Tr stands for trace and Ar({ri}) 

are Lagrange multipliers which ensure the normalization of pedcTi}, {ri}), according to 
Eq.( p.ll]) . Sr is the entropy for the disorder variable r^, for which we also use the CVM 
estimate: 

S-r = Tr^(p^,6lnpr,6) - 3Tr^(p^,2lnpr,2) +2Tr^(p^,ilnp^,i), (3.14) 

where pr,2 and pr,i are partial (symmetrized) traces of pr,6- With the above definitions, 
our free energy can be regarded as a function of pg only and taking the derivative with 
respect to a generic element of pe we find the stationarity conditions 

K . h . 

PeiWi}^ {^i}) = exp[-A^({ri}) + y ^^'^i+i'^^^i+i + 3 5Z 

1=1,6 i=l,6 
X b2(c^l, CTe, n, T6)p2(cri, (72, Ti, r2)p2(cr3, (72, Ta, T2) (3.15) 
P2(cr3, Cr4, T-3, r4)p2((75, (^4, T5, r4)p2(a5, ae, Ts, Te)]^/^ 
X [pi^((Tl)piB((T2)Pl^((73)piB(cr4)pi^((T5)piB((76)]""^/^, 

with the convention (T7 = (Ti, T7 = ti. 



^ The cumulant expansion of the entropy can be understood as follows. We write the truncated 
entropy as 

S= 5h- ^ cSp+ St 

hexagons H pairs P triangles T 

where we must first subtract the contribution of pairs of neighboring triangles and re-add that 
of single triangles to avoid overcounting. Noting that the numbers Nh, Np, Nta and Ntb 
of respectively hexagons, pairs and triangles of the sublattices A and B satisfy Np/Nh = 3, 
Nta/Nh = Ntb/Nh = 1, this leads to the entropy per hexagon appearing in ( ^.13| ). 
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One can solve this set of equations with the definitions of Eqs. (|3.12 ) and the nor 



mahzation constraint by the so-caUed natural iteration method [^. Starting from 



some assumption on p2 and iterating the above equation, pg converges to a solution of 
( p.l5|) which is moreover a local minimum of the approximate free energy (|3.13|) . To find a 



global minimum, it is in general necessary to start the iteration with different sets of initial 
conditions on p2 appropriately chosen to reach the different expected phases. At each step, 
the normalization condition ( 3.11j ) is recovered by adjusting the Lagrange multiplier for 



each realization of the disorder on the hexagon. Before going to the analysis of the model 
with disorder, in the next subsection we will revisit the pure case [|16[. One reason is that 
we need to fix the probability distribution for the disorder variables PriiTi})- 

3.2. Pure case and explicit form for pj. 

In this section, we reconsider the pure case for /i = in detail within the CVM 
approximation as a particular trivial realization of disorder where all Kij are constant and 
equal to K [0. We only need to consider the "pure" 6-point functions p^i^cri), which 
does not depend on the variable any longer. We can easily recognize that the 2x11 = 
22 elements of p^, which correspond to the weights for each state in Fig.l, are not all 
independent since some of the states are related by simple symmetries and should thus 
have the same weight. This of course assumes that the corresponding symmetries are not 
spontaneously broken. Hereafter we only consider the system for /i = and consider three 
types of solutions corresponding to the three different symmetries of the spin system in 
the phase diagram of Fig. 2 [p!6| : 

(1) Disordered folded phase: we do not allow for any spontaneous symmetry breaking in 
the system. Using rotational symmetry and the symmetry under reversal of all spins, 
we end up with only 4 independent weights .^0,1,2,3 corresponding to vertices with 
respectively 0,2,4 or 6 surrounding folds. 

(2) Ferromagnetic phase: we allow for a spontaneous ferromagnetic symmetry breaking 
(M 7^ 0). Then the two vertices with no fold have different weight Zq and Zq according 
to their ±6 magnetization. The other vertices are neutral in this respect and we end 
up with 5 different weights. 

(3) Antiferromagnetic phase: we allow for a spontaneous antiferromagnetic symmetry 
breaking (Mst 7^ 0). Then all weights Zi have to be doubled into (Z^, Zi) except for 
the vertex with no fold [i = 0) which is neutral in the staggered magnetization. We 
end up with 7 different weights. 
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Spin Conf. 


(1) Dis. 


(2) Ferro. 


(3) A.Ferro. 


Deg. 


+ + + + ++ 


Zo 


Zq 


Zq 


1 




Zo 


^ 


Zq 


1 


+ + H 


Zi 


Zi 


Zi 


3 


h ++ 


Zi 


Zx 


'zl 


3 


+ h+- 


Zi 


Z2 


Zi 


6 


- + H + 


Z2 


Z2 




6 


+ - + -+- 


^3 


Zz 


Zz 


1 


- + - + -+ 


Zz 


Zz 


'Zz 


1 



Table I: Independent hexagon spin configurations. 
The corresponding elements of pe and the degeneracies are indicated. 



The spin configurations, their degeneracies and the notations for their weights are 
summarized in Table I. Of course, the case (1) can be recovered from either case (2) or (3) 
as a particular realization with no spontaneous symmetry breaking (i.e Zi = Zi for all i). 
Also, we assume that the two (ferromagnetic and antiferromagnetic) symmetries cannot 
be broken simultaneously. We thus need to study only the cases (2) and (3) above to get 
the complete phase diagram of the system (here at h = 0). 

At first we consider the ferromagnetic case (2). In this case, the stationarity condition 
reduces to the following nonlinear equations between the weights ^0,^0,^1,^2,^3. 

Zo = exp(-A + 3i^)(y++)V(2/+)' 
Z^ = exp(-A + 3K)(y__)V(y-)' 

Zi = exp(-A + K)(y++)(y__)(y+_)/(y+)(y_) (3.16) 
Zi = exp(-A - K)iy++)'/\y__Y/\y+_f/iy+)iy_) 
Z3 = exp(-A - 3K){y+.f/{y+){y.), 

involving a single Lagrange multiplier A. Here, ,y and j/_ are two- and 

one-point functions pi{criai) and pi(cri) (there is no difference here between sublattices A 
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and B), which are defined as foUows, 



y++ = P2(++) = Zo + 2Zi + 2Z2 

y+- = P2{+-) = P2{-+) =Zi + 4Z2 + Z3 

= p2{—) = Z^+2Z, + 2Z2 (3.17) 

y+ = pi(+) = Zo + 3Zi + 6Z2 + Z3 

y- = Pi(-) = Zo + 3Zi + 6Z2 + Z3. 



The above equations imply the foUowing simple relations: 



Zi Z2 Z3 y+ 

Introducing the two reduced variables 

, , y^^. (3.19) 

we can express all the weights in terms of x, y and the normalization factor wq = {ZqZq)^/'^. 

rj +1/2 Ty- -1/2 

Zo = y^ ' Wo, Zo = y ' wq 

(3.20) 

Zi — XWq, Z2—X Wo, Z3—X Wo- 

The above equations ( |3.16|) reduce to the following nonlinear equations for the reduced 
variables, 

_, y + 2y^/^x + x^) .3 A + y^'^{?>x + + x^) . 2 

y 4+2yl/2(x + x2)^ ^y + y^/2{^X + Qx^ + X^y ' ^ ' ' 

^ ^ y{x + Ax'^ + x^) 

w(l + 22/V2(a; + a;2))V2(2^ + 2yV2(a; + a;2))V2- ^ ' > 

where u = exp{2K). The parameter y measures the spontaneous ferromagnetic symmetry 
breaking while the parameter x measures the fugacity per folded bond. 

We can easily see that Eq.( p.21|) has two obvious solutions: a solution y = I and 
X arbitrary and a solution x = and y arbitrary. The latter solution is also a solution 
of Eq.( |3.22]) . It means that each vertex of the membrane can be only in one of the two 
configurations without fold Zq or Zq. The solution cannot determine the ratio y = Zq/Zo, 
i.e the proportion of each state. However, with only these two vertices at hand, no fold 
can be ever created and the only possible global states for the lattice are the state with 
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all spins up (M = 1) and that with all spins down (M = —1). The above solution simply 
describes an arbitrary superposition of these two (symmetric) pure flat states. The fact 
that the membrane is indeed frozen in a pure completely flat state is further confirmed 
by computing the entropy which is found to be exactly zero, and by computing the free 
energy, which is found to be / = —3K per hexagon, as expected (there are 3 bonds per 
hexagon) . 

The first solution with y = 1 means that the spontaneous symmetry breaking does 
not occur and that the membrane is in the disordered folded state (Zq = Zq). The value 
of X is then fixed by Eq.( |3:2^) p 



u 



X 



+ V(3-t^- 



{2u-l] 



(3.23) 



which has a solution for K < ln((l + v 13)/2)/2. Comparing the corresponding free energy 
to that of the pure fiat state, we get a first order transition from disordered folded to purely 
fiat at Kc ~ 0.1013 |T| . 

We also looked numerically for another non-trivial solution with spontaneous symme- 
try breaking {y ^ 1) and intermediate magnetization {x ^ 0) but did not find any. We 
conclude that there is no possible fiat phase with < |M| < 1 and the above three phases 
(M = ±1 or 0) are the only stable ones for positive K. 




2 
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Fig. 4: Probability distribution for each disorder configuration. We also 
show their degeneracies 



As has been discussed previously, the above analysis is also instrumental for the es- 
timation of the probability distribution /9T,6({^i}) of the disorder variables r. We can use 
indeed for p^- g the distribution pe above at -ftT = /i = 0, characterized by y = 1 and 
a: = 2 [|1^. In other words, if we define -Pq, 1,2,3 as the weights Pt-,6 for the local realiza- 
tions of disorder with 0,2,4 or 6 creases around the vertex, we learn that the ratios Pi/Pq, 
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P2/P1 and P3/P2 must all be identical and equal to 2. Their values are then fixed by the 
normalization: 

2Po + 6P1 + I2P2 + 2P3 = 1 (3.24) 

leading to [ 1T6| |: 

— 12 4 8 

= 78'^^ = 78'^^ = 78'^^ = 78- ^^'^^^ 
More generally, we can parametrize the distribution of P with one parameter a equal 
to the ratios a = Pi/Pq = P2/P1 = P?,/P2- Beside the natural value a = 2 above, the 
limiting case a = describes a membrane with no crease and Kij = K everywhere, while 
a = 00 describes a membrane with creases everywhere and Kij = —K on each bond. 

Now we discuss the compactly ordered (antiferromagnetic) folded phase (3). There 
are 7 independent weights for h = (see Table I) and in this case it is convenient to use 
staggered variables 77^ = (— l)*~^ai with (—1)*"^ = 1 on triangles belonging to the sub- 
lattice A and (—1)*"^ = —1 on triangles belonging to the sublatice B. The corresponding 
two-point function is simply 

Pv,2(Vi, m) = P2(?7i, -m)- (3.26) 

About the one-point function, we have the symmetry piA(a) = Pib(— a) in the antiferro- 
magnetic phase, leading to only one (A/B independent) one-point function for 77: 

PnAv) = PiAiv) = Pib{-v) (3.27) 

As before, the solution of the non-linear stationarity equations for the 7 weights can 
be parametrized as: 

r7 rr 1/2 -1/2 

Zq = Wo, Zi = xy ' Wo, Zi = xy ' wo 

rr 2 1/2 IT' 2 —1/2 ry 3 3/2 17~ 3 -3/2 

Z2 = X y I Wo-, Z2=x y ' wo Z^, = x^y"' wq, Z3 = x-'y wq. 
with two reduced variables x and y solutions of. 



(3.28) 



where 



y = ^(^)^/^ x = u-'y^^. (3.29) 

y— y+ y+- 

y++ = Pr,,2{++) = Zi + 3Z2 + Z;+ Z3 

y+- = p^,2(H — ) = Pr,2i — \-) = Zq + Zi + Zi + Z2 + Z2 

y__ = p^,2( ) = Z; + 3Z;+Z2 + Z^ (3.30) 

y+ = pr,,i{+) = Zo + 2Zi + Z; + 4Z2 + 2Z;+Zs 

l/_ = p^,i(-) = Zo + 2Zr + + 4^ + 2Z2 + Z^. 
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In the equations (|3.29| ), the global normalization wq drops out, so the equations can be 
solved for x and y as functions of u = exp{2K). Again the variable x measures the 
fugacity for each fold and y measures the antiferromagnetic spontaneous symmetry break- 
ing. Solving the above equations numerically by iteration, we find a continuous transition 
from a disordered folded state (y = 1) to a compactly folded ordered phase (y 7^ 1) at 
Kst = —0.2838. The value of Kgt can be found simply by linearizing the equations ( |3.29| ) 
by writing y = 1 + e. This fixes the value of x to be the real solution Xgt of 



x2-21x2-12x-4 = (3.31) 



that is 



and u to be 



Xst = 7+ (387 + 2^223)^/^ + 53/(387 + 2^223)^/^ (3.32) 



1 + 4xst + xj^ 
= l + 2x., + 2xl- 



4. Results for the fully disordered system 

In this section we analyze the fully disordered case (Eq.( p.9| )) within the CVM ap- 
proximation. In subsection 4.1, we study the system for h = for several values of the 
parameter a for the disorder weights. Next we fix a = 2 and proceed to the general (-ftT, h) 
case in subsection 4.2. We obtain the {K, /i)-phase diagram by use of the natural iteration 
method. 



4.1. Analysis with reduced elements of pq (h = case) 



As in the previous section, we will study the fully disordered system for 
/i = by reducing the number of elements of the 6-point density matrix. In the case 
with disorder, the symmetries of the elements of the 6-point function pi5{{(Ti}, {ri}) de- 
pend also on the symmetry of the disorder (r- variables) , in addition to the symmetries of 
the spin variable {a} itself. 

In the previous section, we have studied the pure model with an antiferromagnetic 
spontaneous symmetry breaking for K < 0. There we have used the staggered variables as 
r]i = {—iy~^ai. If we regard the pure system for K < as a trivial disordered system with 
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xconfig. T^config. x config. T^config. 




Ticonfig. 



Fig. 5: The 38 independent local fold environments for each vertex for the 
4 different local realizations of disorder (0,2, 4, or 6 creases). Each disorder 
configuration is shown at the left hand side of each group. To their right, 
we show the spin configurations, the weights and the degeneracies. We have 
represented by thick lines the domain walls for the gauged spin variable ?7i = 
(TjTj. A subscript i,j indicates a configuration with 2i creases and 2j folds. 
The superscript ± indicates a ± contribution to F. In a F = phase, equating 
the + and — weights leaves us with 22 independent weights. 

Pt,6{^i —1, 1, —1, 1, ~1) = 1, i-e. Ti is fixed to (—1)*"^, and K > 0, the above staggered 
variables can be written as 

r]i = GiTi. (4.1) 

The motivation for introducing the staggered variables is that in these variables the anti- 
ferromagnetic order parameter Mgt is simply written as 

= (4-2) 
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and we don't have to differentiate between the A- and B-sublattices. That is, the staggered 
variables {77^} are more natural than the original variables {a^} when one discusses the 
antiferromagnetic symmetric case. 

In the fully disordered system, we are mainly interested in the spontaneous symmetry 
breaking at iiT > of the following "frozen" order parameter! 



This order parameter judges whether or not the membrane is trapped in the randomly 
oriented phase, characterized by the disorder variables {r^}. As in the antiferromagnetic 
case, it is natural to use the following "gauged" variables. 



Using these new gauged variables, we classify the elements of 6-point functions by the 
symmetries of both the spin configuration and the disorder configuration. Hereafter we 
only consider the system for h = and allow for two types of solutions which correspond 
to whether the frozen order exists [F 7^ 0) or not (F = 0). 

(1) Disordered folded phase: we do not allow for any spontaneous symmetry breaking 
in the system (M, Mst and F = 0). Each of the 4 elementary types of disorder 
configuration (with 0,2,4 or 6 creases) leaves us with a certain number of symmetries, 
including that under reversal of all spins. We use these symmetries on the 77 variables 
to reduce the number of weights. We end up with only 22 independent weights in this 
case. 

(2) Frozen phase: we allow for a spontaneous symmetry breaking of the frozen order 
parameter [F 7^ 0). Then all weights have to be doubled except for those vertices 
which are neutral in the gauged magnetization (YliVi = 0)- We end up with 38 
different weights in this case. 

In Fig. 5, we have summarized the results of this symmetry analysis in the case (2). 
Note that case (1) can always be seen as a particular case of case (2) with extra symmetries. 
On the left hand side of each group, we show the disorder configuration {r^}. To its right, 

^ The reader might wonder whether the r ^ — r symmetry could lead to a zero quenched av- 
erage of the order parameter. However, this symmetry implies only that whenever /Omin({o"}, {t}) 
is a solution of the variational equations, Pmin({<7}, {— r}) is also a solution, but does not imply 

that Pmin({cr}, {r}) = Pmin({0-}, {-t}). 




1 



(4.3) 



(4.4) 
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^0,0 ^0,0 %,i 




X,0 X,! ^1,1 

X config. r| config. 

Fig. 6: Definitions of two-point functions. Gauged spin variables rji = aiTi 
are used. On tlie left hand side, we show the disorder configuration r. 

we present the spin configurations {rji} which are independent from each other. We also 
indicate the notations for their weights and their degeneracies. The two indices i,j in Zij 
indicate a configuration with 2i creases and 2j folds in the a variable. 



In terms of these elements of the 6 point function, we define the two-point functions 
as follows, 

yo,o ~ -^0^0 + 2Zo,i + 2Zo,2 + 2Zi^Q + 2Zf^-^ + 2Zf -^ + 2Z^2 
+ '2,Z^2 + 2^2,0 + 2.^^! + 2Z^]^ + 2Z^2 + 2-^2^2' 

2/0,1 = -^0,1 + 4Zo,2 + -^0,3 + Z^^i + -^1,1 + Z^^2 + ^1,2 
+ 3Z^2 + 3^1,2 + ^1^3 + ^1,3 + ^2^1 
+ -^2,1 + -^2,2 + -^2,2 + 2-^^2 + 2-^2^2 

+ -^2^2 + -^2,2 + -^2^3 + -^2,35 (4-5) 

yti = + zt,2 + ZI2 + ^zt2 + zt^ + Z2±i + zl^ 

+ 4Z2±2 + 2Z^2 + 44t2 + 2% + 4412 

+ 3-^2^3 + -^2^3 + 2Z^i + Z^3 + 3Z^2 
+ Z^2 + -^3^1 ■ 

The upper right index of y means that the gauged spin configuration {r\i\ have the cor- 
responding positive (or negative) contribution to F. Again, the first lower right index 
indicate whether there is a crease line (1) or not (0). The second index means that there 
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is a fold (1) or not (0) in the original a variables. Based on these two-point functions, we 
also introduce the following one point functions, 

^Vo,Q + yo,i+yi,o + yti- (4.6) 

In terms of these functions, we write down the stationarity conditions. For example, 
let us show those for Z^q, Zo,i, Zo,2 and Zq^z- 

Zo+o = exp(-Ao + m){yl^f/{y+)\ 
Z-^ = expi-Xo+3K){y-,)y{y-f, 

Zo,i = exp(-Ao + K)iy+,)iy-,)iyo,i)/iy+)iy-), (4.7) 

Zo,2 = exp(-Ao - i^)(yoto)'/'(yo-o)'/'(yo,i)V(y+)(y-), 

Zo,3 = exp(-Ao - SK){yo,if/{y+)iy-). 
Here, due to the above symmetries, we need only 4 Lagrange multipliers Aq, 1,2,3, one for 
each of the 4 elementary types of disorder in Fig.4. These Lagrange multipliers are of 
course determined by the normalization conditions, like 

Zo,o + Zo,o + 6^0,1 + 12Zo,2 + 2Zo,3 = Po- (4.8) 

As in the pure case, we introduce reduced variables. Here we need 4 ratios x, y, s and t 
defined as. 



-1 yo,i -1 yi,o 

x = u , , _ y = u 



(yo>o":o)^/'' (ytiyi,!)'/'' 
/yJoxi/2.y-xi/3 t= (^M\V2/y-xi/3 

S/o"o y+ ' yi,i y+ 



(4.9) 



(4.10) 



and for convenience, we also introduce the following averaged weights, 

wo = {Zl^Zl^fl\ w, = {Zl^Zl^fl\ 

W2 = (^2+2^2:2)'/', m = {zlzZi^Y^\ 

The stationarity conditions are then reduced to the following simple form, 

= 'S'^^^'o, ^0,1 = a:wo, ^0,2 = x^wo, ^0,3 = x^wq, 

rv 3 r7± 2.±1 ^± .±1 ^± j.±3 

^3,0 = y W3, Z3 1 ^y t ^3, Z3 2 = «;3, ^3,3 = ^ ^3, 
■^1,0 = y^x, -^1,1 = s t tui, 1 = xys wi, 

= xs^'wu Z^,2 = x^/V^Hst)^^/''wi, (4.11) 

r7± 2.±1 ^ 2 ^± ±1 

3 = X t Wi, Z2,o = y W2, Z2 1 = ys W2, 

= x^/2y^/2(st)±^/^i(;2, ^^2 = S^^t^^«'2, ^^2 = xyt^^W2, 
Z^^^ = xyt^^W2, = a;^/2yi/2(st)±V2^2, ^2^3 = a;t±i«;2. 
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Each weight j is given as the product of Wi by a simple function of x, y, s and t. The 
rules for the x and y variables are simple: The disorder configuration splits the bonds into 
those which support a crease and those which do not. On the bonds with no crease, we 
assign a factor a/x if the bonds has a fold and 1 otherwise. On the bonds with a crease, we 
assign a factor y/y if the bond has no fold and 1 otherwise. In both cases, the non trivial 
factor is assigned if the gauged variable changes sign when crossing the bond. About the 
factors of s and t, the rules are more subtle. Still s and t both measure the symmetry 
breaking of the frozen parameter. 

The weights wo,i,2,3 can be expressed as functions of x, y, s and t thanks to the nor- 
malization conditions as follows, 

Po = wo[s^ + + 6x + 12x^ + 2x^], 

Pi = wi [2y + sH + s~'^t + 2xyz + 2xyz~^ + 2xy + 2xy~^ 

+ 4x3/2^1/2,1/2^1/2 ^ 4^3/2^1/2^-1/2^-1/2 ^ ^2^ ^ ^2^-1]^ 

P2 = W2[2y^ + ys + ys-' + 2x^2^^/2^1/2^1/2 ^ 2^1/2^3/2^-1/2^-1/2 ^^ -^^^ 

+ st^ + + ?,xyt + ?,xyt-^ + 2x^/22/i/2si/2^i/2 

^2x'/V^s-'/h-^'^+xt + xt-'] 
P3 = wslt"^ + + Qy{t + r^) + 32/2(t + + 2y\ 

In Eq. (|4.9|) , the right hand side of each equation is thus a function of x, j/, s and t only. 
These equations can be simply solved numerically by iteration. 

Hereafter we show the results of this numerical analysis. We first fix a = 2 for the 
weights -Po,i,2,3- We start the iteration for K = Q with a fully symmetric solution which 
corresponds to a disordered phase. We proceed to the iteration until the required precision 
is reached. We then increase K by dK and restart the iteration. For this next value of 
K, we start the iteration from the solution of the iteration for the previous value of K. 
This procedure allows to follow the continuous evolution with K oi a, given local minimum 
of the free energy. We increase K from to 0.2 and then decrease it back to 0. In this 
way, if the system has a first order transition with two local minima of the free energy 
in competition, the method will show a hysteresis. Note that this iteration procedure 
is slightly different from the natural iteration method that we shall use in section 4.2, 
where we search a solution for the nonlinear stationarity equations from different initial 
assumptions corresponding to the different possible symmetries. 
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Fig. 7: Frozen order parameter F versus bending rigidity K. F changes from 
F = to F 0. The system shows a hysteresis with two separate jumps 
for two values of K on each side of Kp- The position of Kp is determined 
precisely by comparing the free energies of both phases, as shown in the next 
figure. 

In Fig. 7, we show the behavior of the frozen order parameter F as a function of the 
bending rigidity K. For small K{< Kp), F is clearly zero and the system is in a disordered 
folded phase. At Kp, the system shows a first order transition from this disordered folded 
phase to a frozen phase F ^ 0. The value of |F| is strictly less than 1, thus the system 
is only partially frozen. When the iteration is performed with first increasing K and then 
decreasing it back to zero, we see a clear hysteresis with two jumps on both sides of Kp. 
The value of the transition point Kp can be fixed precisely by comparing the value of the 
free energies for both phases. This is shown in Fig. 8. The hysteresis allows us to see clearly 
the crossing of the two free energy lines corresponding to both phases F = and F ^ 0. 
Indeed the system stays for some time after the transition point in the wrong metastable 
state. As shown in the inset of the figure, the transition occurs at Kp ~ 0.166(1). 

Fig. 9 shows the behavior of the two-point function < '?7i?72 > = —internal energy /K 
versus K. It also shows a clear evidence of first order transition with a hysteresis in 
the results of the iteration procedure. At K = 0, the disorder variables Tj and spin 
variables decouple. The value of the function is then < 771772 > = T1T2 < cricr2 >atK=o= 
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-0.65 ^ ,x ' 

0.00 0.05 0.10 0.15 Kp 0.20 

K ^ 

Fig. 8: Free energy / per hexagon versus bending rigidity K. The thin 
straight hne corresponds to a completely frozen phase F = 1. This phase 
is never stable. The two thick lines correspond to disordered folded phase 
F = (for small K) and to the frozen phase < |F| < 1 (for larger K). The 
two lines cross at the transition point Kp. As shown in the inset, we find 
Kp ~ 0.166(1). 
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Fig. 9: Gauged two point function < ?7i?72 > = —internal energy/iiT versus 
K. We see here also a clear evidence of first order transition with a hysteresis. 



(—1/3) X (—1/3) = 1/9, as found here (the value —1/3 is easily obtained from the analysis 
of the pure case at iiT = of section 3.2 ill 
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Fig. 10: Frozen order parameter F versus K for several values of a (we only 
show here the case -F > 0). The jump in the order parameter becomes smaller 
as a becomes larger. The continuous character of the transition is recovered 
at a = oo (Pure antiferromagnetic system. For intermediate a, the iteration 
gives rise to a hysteresis.) 

Next, we studied the system for several values of the parameter a for the disorder 
weights. As has been explained previously, a = means that there is no crease in the 
system and a = oo corresponds to the pure antiferromagnetic system. From the previous 
analysis, we know that the system shows a first order transition for both a = and a = 2 
above. At a = (where F=M) the transition is from F = M = OtoF = M = ±l |T^ . 
At a = 2, the discontinuity is smaller with |F| < 1 in the frozen phase. We also know that 
the transition becomes continuous and second order at a = oo [16|. The transition point is 
aX K = 0.284 as obtained before. From the results in Fig. 10, we see that the discontinuity 
of the transition becomes smaller as we increase the parameter a. The continuity of the 
transition seems to be recovered only at a = oo, although it is difficult to determine 
whether the transition is of first order or of second order when the discontinuity becomes 
too small. 



4-.2. (K,h) -Phase diagram 



Let us now turn to the analysis of the whole phase diagram in the {K, /i)-plane. It is of 
course symmetric with respect to the /i = axis. It is shown in Fig. 11 for /i > as obtained 



from the CVM stationarity equations solved by the natural iteration method [21 1. Here 
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Fig. 11: Phase diagram in the {K,h) plane for the fully disordered system 
(Eq. (|2.9|) ). First order lines separate the four phases: (1) Disordered folded 
phase with M ~ 0, F = 0, (2) Completely Flat Phase with M = 1 and F = 0, 
(3) Flat Phase with < |M| < 1 and F = and (4) Frozen phase with M 
and |F| > 0. 

we have set a = 2 again and the results were obtained with a set of 22 x 22 independent 
weights, i.e. without making any assumption on the symmetries of the different phases, 
except in the initial conditions of the iteration. At sufficiently large values of h and 
for h > 2K, the completely flat phase with M = 1 is stable with respect to both the 
disorder and the thermal fluctuations. At sufficiently large values of K and for h < 2K, 
the system is in a frozen phase with M = 0, Mgt = and F ~ 1. There is a first 
order transition line between these two phases, which is roughly given hy h = 2K. The 
position of the line can be obtained by requiring that the free energy of the completely 
flat phase and that of the frozen phase take the same value, i.e. by solving the equation 
/c.flat = /frozen- About /c.flat, the bending energy contribution per hexagon is estimated 
as — 3i^r< rir2(Jicr2 > = — 3i^r(rir2) = —?)K x (—1/3) = K, while the entropy vanishes in 
the absence of local excitations. We thus get the exact free energy /c.flat = ~2h + K. The 
estimation of /frozen is more difficult and we simply assume that /frozen — — 3-ftr as in a 
completely frozen phase, because the frozen order parameter F is almost saturated to 1. 
From these estimations, we obtain the transition line h ~ 2K, which is what we indeed 
observe. 

For smaller K and h, the system is in a disordered folded phase with F = and 
M ~ 0. As for the pure system, M does not vanish exactly for h > but still remains 
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very small. This might be an artifact of the CVM approximation [^. At i^T = 0, the spin 
variables a are decoupled from the disorder variables r and the fully disordered system is 



the same as the pure system. It has a first order transition point at /i ~ 0.184 [jT5|. For 
/i = 0, the above results confirm those of previous section with a transition aX K — Kp- 



Disorder 
Config. 







Fig. 12: Lowest energy spin configurations {0"^} for each type of disorder 
configuration. Folds are indicated by thick lines and creases by dashed lines. 
On the left hand side, we show the disorder configurations, in the center the 
corresponding lowest energy state, which violates the local folding constraint 
for the disorders with 1 and 2 creases, and on the right hand side the lowest 
energy states which preserve the constraint, together with their degeneracy. 



We have also studied the case of negative although it is not very physical. Still, 
it presents interesting features in view of our further study of other types of disorder in 
Section 5. In the usual Mattis Model without constraint and for /i = 0, the spins develop at 
large enough negative K an "antiferromagnetic"-like order in the gauged variable rji with 
a ground state rji = (—1)*"^ in the limit K = —00. Here such order cannot be reached 
in general due to the constraint on cr, hence on 77. If the disorder has or 3 creases, 
then the ground state rji = (—1)*"^ can be reached and is the unique ground state spin 
configuration. On the other hand, if it has 1 or 2 creases, it cannot be reached and we 
are led to several lowest energy spin configurations, as shown in Fig. 12. At K = —00 
and h = 0, the actual ground state will thus be degenerate with frustrations in the system 
which might prevent the emergence of a true "frozen antiferromagnetic" order. The system 
is thus always disordered in contrast with the pure case where an antiferromagnetic order 
had developped. 
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Fig. 13: Gauged two point function < ?7i?72 > = <~dYfia2T2> along the 
K-axis {K <0,h = 0). 

From the analysis of the lower energy states, we can easily compute the two-point 
correlation at K — > —oo, c = < aiTia2T2 > for the gauged variable (i.e. the internal 
energy /K). For a disorder with or 3 creases, we have J2i around v ^ (^i'^if^i+i'^i+i >= — 6 
in the lowest energy state. For a disorder with 1 or 2 creases, we have around v < 
crjTjCTj+iTj+i >= —2 for all the lowest energy states. The averaged value is thus estimated 
as c = < c7iTi(72T2 > = (-6 (Pq + ^'s) " 2 (Pi + P2))/6 = -19/39 ~ -0.48. This is what 
we observe (see fig 13.). The two-point function does not show any discontinuity in all the 
negative K regime, which confirms the absence of transition in this regime (at h = 0). 

Finally, we included in Fig. 11 the results of the CVM analysis for negative K and 
arbitrary h. We can see the emergence of a new partially flat phase with < |M| < 
1. The nature of the transition between this phase and the disordered folded phase is 
unclear, in particular because in the latter phase, M is not exactly zero within the CVM 
approximation. We then see a limiting point (black circle in Fig. 11) below which the 
magnetization does not present any longer a discontinuity between the two phases. This 
also might be an artifact of the CVM, in which case the true transition line should be 
continued to lower values of K (dashed line). 

5. Other Models 

In this section, we complete our study by considering other variants of the disorder. 
As discussed in section. 2, there are several possibilities for the choice of Kij. In order to 
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appreciate the importance of the local folding constraint (|2.1| ) on the disorder face variables 
Ti, we will study the model Hamiltonian ( p.9|) without the local folding constraint on r^. 
Next we will study the Edwards- Anderson model with the local folding constraint on the 
spin variables ai and a bending term Kij given by Kij = Krij with a random variable 
Tij = ±1 on each bond (ij). We refer to the former case as the model (2) and to the latter 
case as the model (3). The difference between models (2) and (3) is simply the possibility 
in the model (3) of having vertices with an odd number of surrounding creases. We also 



refer to the original model ( |2.9|) with the folding constraint on as the model (1). 
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Fig. 14: One of the lowest energy spin configurations {ai} for each type 
of disorder configurations. Folds are shown by thick lines and creases are by 
dashed lines. The degeneracies for the disorder configurations are indicated. 
We also give the value of the corresponding minimal energy. 

Before we discuss the corresponding {K, h) phase diagrams, let us discuss the two-point 
function of the system at /i = and large K. Clearly in this limit, for a fixed disorder 
configuration, the ai variables will tend to minimize the energy —Kijaiaj, i.e. will tend 
to maximize the overlap with the disorder configuration in terms of folded bonds. In 
other words, the system wants to create a fold (ajCTj = —1) whenever a crease exists 
{Kij/K = —1), and no fold otherwise. For an arbitrary environment of Kij around a 
vertex, we can easily find one corresponding lowest energy state for the a variables around 
the vertex. There are in general several such states. In Fig. 14, we have displayed all 
disorder environments together with one of the corresponding lowest energy state. It is 
interesting to notice that the disorder configurations can be arranged in three categories: 
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(i) Those with an even number of creases and which satisfy the folding constraint. The 
lowest energy state is unique and has energy —6K. 

(ii) Those with an even number of creases but which do not satisfy the folding constraint. 
The corresponding minimal energy is —2K in this case. 

(iii) Those with an odd number of creases. The corresponding minimal energy is —4K in 
this case. 

Of course, in the model (2), only vertices of type (i) and (ii) are allowed while in the 
model (3), all vertices can appear. 

At — > oo and h = 0, we can thus estimate the two-point function < 771/72 > = 
< (JiTia2T2 > for the model (2) and < criTi2(T2 > for the model (3) by averaging over 
all disorder environments the corresponding minimal energy. All (allowed) disorder envi- 
ronments are now equiprobablc. Taking into account the appropriate degeneracies under 
rotations, we get for the model (3) 

< (^1X12(72 > = (-6x1-6x1-4x6-2x6-2x6-6x3-4x6 

(5.1) 

-4x2-4xl2 -2x3-2x6-6x6-4x 6)/(6 x 64) = 59/96 - 0.614. 
For the model (2) we find easily: 

< (J1T1CT2T2 > = 18/32 - 0.562. (5.2) 



0.80 




K 



Fig. 15: Gauged two-point function < criricr2T2 > for the model (2) (solid 
line) and < criTi2cr2 > for the model (3) (dashed line) versus K for h = 0. 
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Note that the above calculation assumes that a lowest energy state can be constructed 
globally out of these local lowest energy configurations. This assumption is acceptable 
within the CVM approximation at least. We show on Fig. 15 the two-point function for 
the models (2) and (3) as obtained from the CVM They do not display any discontinuity 
and tend to the values calculated above at large K. 

Note also that in some sense, the model (3) is less frustrated than the model (2) since 
a better overlap with a constrained a configuration can be obtained for those frustrated 
disorder environments with an odd number of creases. 
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Fig. 16: Phase diagram in the [K, h) plane for the model (2). We find three 
phases: (1) a disordered phase with F = g = 0, (2) a completely fiat phase 
with M — 1 and F = q = 0, and (3) a frozen fiat phase with < M < 1, 
q > and F > 0. The solid line represents a first order transition line and 
the dashed line a continuous transition line. 

We have studied the phase diagrams of both models (2) and (3) within the CVM 
approximation by use of the natural iteration method |2^. In order to characterize each 
phase, in addition to the previous order parameters M and F we will also use the spin-glass 
order parameter q defined in Eq.( |2.7| ), generalized to the case where M ^ 0: 

i 

In the model (2), both F and q measure the frozen character of the phase. They are 
expected to be zero (resp. non-zero) simultaneously for our choice of disorder distribution. 
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Fig. 17: Phase diagram in the {K,h) plane for the model (3). Two phases 
are separated by a first order transition line: (1) a frozen fiat phase with 
< M < 1 and q > 0) and (2) a completely flat phase with M = 1 and q = 0. 
The disordered folded phase with M = q = is recovered only at h = 0. 

In the model (3) with no face variables, F is not defined any longer and we will use q 
as a measure of the frozen character of the phase. The phase diagram for the model (2) is 
shown in Fig. 16 and that for the model (3) in Fig. 17. They are of course symmetric with 
respect to the K = axis but also with respect to the h = axis. This is because the 
transformation — > (— l)*~^Tj in the model (2) and — > — in the model (3) exchange 
equally probable disorder environments and simply change K into —K. We only display 
the phase diagrams for K > and h > 0. 

At K = 0, both systems are identical to the pure system and undergo a first order 
transition at /i ~ 0.184 to a completely flat phase. This completely flat phase with |M| = 1 
and = is stable for all K at large enough h above a line which is almost the same for 
the two models up to the tricritical point of Fig. 16. At h = 0, both system remain in a 
disordered folded phase with M = q = for any value of K. This is different from the 
model (1) where we had a transition to a frozen phase at K = Kp- The absence of a frozen 
phase at /i = is again due to the presence of frustration in the system leading to several 
competing lowest energy states. 

For fixed (large enough) K and increasing h, the models (2) and (3) display different 
behaviors. As far as q is concerned, the model (3) develops a non-zero value of q for any 
h > 0. On the other hand, q remains zero in the model (2) until a critical value of h is 
reached where a continuous transition to q ^ occurs. For both models, the "frozen" 
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phase with q ^ also has < |M| < 1, and is thus partiaUy flat. As to the q = phase 
of the model (2), we see a non-zero value of M which is indeed non negligible close to the 
continuous transition line. Still we cannot exclude that this could be again an artifact of 
the CVM approximation. Indeed, by continuity from the K = line, we would rather 
expect M = everywhere in this phase. This issue is thus not fully solved. Finally, the 
absence of the q = phase in the model (3) (except for h = 0) might also be interpreted 
again as an indication of a weaker frustration as compared to the model (2). 

6. Discussion and Concluding Remarks 

In this paper we have studied the folding of the triangular lattice in the presence of a 
quenched random bending rigidity Kij = ±K and a magnetic field h. We have considered 
three types of quenched randomness (1) Kij = KriTj with face random variables ± 1 
subject to the folding constraint ( |2.1| ); (2) Kij = KriTj without the folding constraint 
on the Ti's; (3) Kij = Krij with a bond random variable Tij = ±1. In case (1), the 
folding constraint on the disorder variables was introduced to describe a particular type 
of "physical" disorder supposed to mimic that induced in a randomly crumpled surface, 
here in the context of the folding of the triangular lattice. Applying the cluster variation 
method generalized to random systems, we have studied the phase diagrams of the three 
models (1),(2) and (3) and their phase transitions. The phase diagrams for each case are 
depicted in figs. 11, 16 and 17 respectively. The most important difference between the 
model (1) and the models (2) and (3) is that, in the absence of magnetic field, a frozen 
phase is found only in the model (1), for large enough K. In this phase, the configuration 
of the triangular lattice is trapped in the randomly oriented phase characterized by the 
configuration of the disorder variables {r^}. The models (2) and (3) do not present such 
frozen order at /i = 0. Indeed, these models, where the quenched randomness is not 
constrained, have strong frustrations because the constrained spins describing the normals 
to folded configuration fail to be in the "virtual" lowest energy ground state dictated by 
the unconstrained disorder, even if the coupling constant K becomes large. For /i > 0, a 
frozen phase is recovered in the models (2) and (3). We find several first order or continuous 
transition lines between the frozen phase and a completely fiat phase or a disordered folded 
phase. 
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At last we make one comment about previous studies on another spin model describ- 



ing a polymerized membrane with quenched random spontaneous curvature [p7| , p8[| , with 
Hamiltonian 

H = — KHi ■ fij — Dij ■ [Hi X fij). (6.1) 

Here n denotes the normal vector to the membrane embedded in 3-D-space. The first term 
is a bending rigidity term and K is the bending rigidity modulus. The second term is a 
random spontaneous curvature term with a Gaussian probability distribution for D^j with 
variance 

< Dl >= r^. (6.2) 

In particular, it does not satisfy "physical" constraints of a spontaneous curvature 
which would have been induced by crumpling. Within a mean field study, it was concluded 
in p7|J2^ ] that the model has a wrinkled phase in the {K, r)-plane with nonzero spin glass 
order parameter q ^ 0. This is to be contrasted with our results where the existence of 
such a phase was crucially requiring the "physical" constraint on the disorder variable. 
However it is not yet clear whether our conclusions for a quenched random rigidity are 
applicable to the quenched random spontaneous curvature case. To study the folding of 
the triangular lattice with random spontaneous curvature, we would need to go to a three- 
dimensional embedding space. One possibility is to introduce disorder in the 96-vertex 
model of ref.|T^. This is left for future study. 
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